No classical limit of quantum decay for broad states 

N. G. Kelkar and M. Nowakowski 
Departamento de Fisica, Universidad de los Andes, 
Cm. IE N0.I8A-IO, Santafe de Bogota, Colombia 

Abstract 

Though the classical treatment of spontaneous decay leads to an exponential decay law, it is well 
known that this is an approximation of the quantum mechanical result which is a non-exponential 
at very small and large times for narrow states. The non exponential nature at large times is 
however hard to establish from experiments. A method to recover the time evolution of unstable 
states from a parametrization of the amplitude fitted to data is presented. We apply the method 
to a realistic example of a very broad state, the a meson and reveal that an exponential decay is 
not a valid approximation at any time for this state. This example derived from experiment, shows 
the unique nature of broad resonances. 

PACS numbers: 03.65.Ta, 03.65.Xp, 13.25.Jx 



1 



I. INTRODUCTION 



It is well known by now that the exponential nature of the decay law which appears 
practically in every field of physics and which follows from classical physics is an approxi- 
mation [l( an d deviations from the exponential are expected at extremely short and large 
times. Apart from the power law behaviour at small times, the quantum mechanical sur- 
vival probability of an unstable state typically displays three regions: an exponential decay 
law followed by an oscillatory behaviour corresponding to the transition region and finally 
a power law behaviour (the non-exponential tail). Over the years there have been several 
unsuccessful attempts in particle and nuclear physics [2J to verify the non-exponential (NE) 
tail experimentally. The failure of such experiments is due to the fact that the critical time 
for the transition from the exponential to the NE depends on the width of the state. For 
narrow states (i.e. small width and hence long lived) the critical time is large (up to several 
years for heavy nuclei) and the exponentially decaying sample would physically diminish to 
an unmeasurable amount. For broad states the critical time is small. However, the dominant 
decay law exp(— Ft) at small times (following the extremely short time region) reduces the 
sample rapidly due to the large value of T in the exponential. In essence, one could say that 
nature conspires to hide the NE tail. 

In the present work, we investigate the survival probabilities of broad states which display 
a peculiar behaviour contrary to the standard picture mentioned above. The scattering 
data in reactions where an unstable intermediate state is formed, can be used to evaluate 
the surviva! probability of the —ate state at all thnes & Using this m ethod which 
leads to an indirect measurement of the NE decay at large times, we show that for a very 
broad state, the decay law is never exponential. The classical limit confirmed so often in the 
laboratory for narrow states, does not exist for broad unstable states. An understanding of 
the time evolution of unstable (sometimes referred to as metastable or resonant) states is 
of fundamental importance for every branch of physics where decaying states appear. As a 
matter of fact, except for the electron and proton every other elementary particle is unstable 
and decays spontaneously. 

Since the main objective of the present work is to present a semi-empirical method to 
evaluate survival probabilities and then discuss the special case of broad resonances, in what 
follows, we shall first introduce the formalism used to evaluate the survival probability. We 
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shall then present a realistic parametrization of pion-pion elastic scattering data and use it 
to demonstrate the result mentioned above for broad resonances. 



II. FORMALISM 

The quantum mechanical survival probability of a decaying state, without any approxi- 
mation is given as, 

Pit) = \A(t)\ 2 = \{^\e~ im \^)\ 2 . (1) 

At very small times, it can be shown to be P(t) ~ 1 — (A^H)t 2 , where, A^H is the 
uncertainty in energy. A direct deviation from the exponential decay law at short times has 
been experimentally verified 4j. This deviation is expected at extremely short times and 
will not be the topic of concern in the present work. The exponential fall of P(t) which 
follows is the behaviour most commonly verified in the laboratory. With the exception of 



the ^ with orgaoic — Q, the larg e thne heh™ has however not oeen 
measured and there exist different theoretical approaches for the evaluation of Pit) at large 

n 

times pj. An interesting discussion on the difficulties and possibilities of measuring the 
non-exponential tail can be found in jf]]. 

A. Fock-Krylov method 

In the present work we use the Fock-Krylov method j^l, [?J which relies on basic results in 
quantum mechanics and is briefly presented below. Given the fact that an unstable state \^>) 
cannot be an eigenstate to the Hamiltonian {H\ty) ^ E\^>) otherwise A(t) = (ty\e~ lHt \ty) = 
e -iEt an j p(t) _ i i m plying that the state never decays) one can expand 1^) as, 



|*) = / dEa(E) \E) , (2) 

JSpect(H) 

where, H\E) = E\E). Using (E'\E) = 5(E' — E), we arrive at the result that 

P (E) = dPT °^ {E) = \a(E)\\ (3) 

is a probability density (and as such positive-definite) {3} to find the states with energy E 
in the resonance. One can now evaluate the survival amplitude, 

/•CO 

-lEt I 7 7-1 „/77T\„— iEt 



A(t)= dE p(E)e- lht = dEp(E)e- lht (4) 

•/Spect(H) JE th 



which turns out to be a Fourier transform of the spectral function p(E). E t h is the sum of 
the masses of the decay products. The general form of p(E) = (Threshold factor) x (Pole) x 
(Form factor), i.e., 

p{E) = (E — E th y x P{E) x F(E) . (5) 

P(E) has a simple pole at zr = En — iTn/2 which leads to the exponential decay law. F(E) 
is a smooth function which should go to zero for large E. Going over to the complex plane 
and performing the integral as described in [3J, the survival amplitude A(t) is given as, 

poo 

A(t) = Res[p(z), z R ] + e~ iEtht / dx P(-ix + E th ) F(-ix + E th ) x^ e~ xt (6) 

Jo 

= A exp {t) + A L {t) 

with A exp (t) and Ai{t) representing the exponential and the remaining part of the amplitude 
respectively. We shall now proceed to evaluate A(t) for a particular choice of p(E) which 
connects it to scattering data. 

B. Density of states 

The connection between scattering data and p(E) as noticed in 0] is briefly repeated 
here for clarity. While calculating the virial coefficients in the equation of an ideal gas, 
Beth and Uhlenbeck js] found that the difference between the density of states (of scattered 
particles) with interaction dni(E)/dE and without dnf ) \E)/dE is, 

dn _ dm(E) _ dnf\E) _ 21 + 1 d5 t (E) 

dE~ dE dE ~~ dE~' ^ 

Si(E) is the scattering phase shift for the I th partial wave in elastic scattering. For an 
intermediate unstable state occurring in the scattering of two particles, this is the density of 
states of the unstable state (or resonance) in terms of the decay products. Thus the spectral 
function p(E) = dPTobq,(E)/dE oc dn/dE can be expressed in terms of the derivative 
of the scattering phase shift (which in turn is related to the scattering amplitude 7] as, 
T\ = [exp(2i5i) — 1]/ 2i). This phase shift derivative was in fact found to be the delay time 
(or phase time delay) in scattering by Wigner [h]] and also used to characterise resonances 



in hadron scattering llj]. This interpretation works well for all /- values except for the 



s-wave (1 = 0), because in this case dS/dE oc (E — Eth) an d we encounter a threshold 
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singularity 12J. The problem can however be resolved by defining rather a dwell time delay 
(which has also been shown to be a density of states 13|) as proposed in 14J. Thus, 



' dn\ , „ , , d5 2Ue(T)-Js 

—— = dwell time delay = 2—— — 

dE J dE s-E? 

\ / new 



(8) 



which is the relativistic version of the expression found in [14(. Here s = E 2 and E t h is the 
sum of the masses of the decay products. With (jSJ) as p(E) one can check that one gets the 
standard threshold behaviour and replacing this p(E) in @ the correct power law as also 
found in Il5l. 



III. THE CASE OF A BROAD STATE: a MESON 



Being equipped with the theoretical framework for evaluating A(t) and hence the survival 
probability P(t) = \A(t)\ 2 , we now proceed to calculate P(t) for a realistic broad unstable 
state. The choice we make is that of the scalar meson a formed in pion-pion (tt it ) elastic 
scattering. The very short lived a has been and is still one of the most controversial prob- 
lems among particle physicists. It was removed from the particle data listing in 1974 and 



reappeared there much later, 
in different physical situations 



t is sometimes claimed that this meson behaves differently 
16j . i.e., displaying different masses and lifetimes. Theoreti- 



17j after 



cally, it can be viewed as a Higgs particle in the context of the linear sigma model 
the spontaneous breaking of chiral symmetry. It can also be looked upon as a low energy 
manifestation of the scale invariance breaking in the strong interaction 18]. Some recent 
discussions on this enigmatic scalar meson can be found in 19]. 



A. Parametrization of the amplitude using tttt scattering data 



To evaluate the density of states for the sigma meson, we use a parametrization of the 
scattering phase shift given in [20] which is obtained from a consistent fit to the production 
and elastic ir it scattering data and includes the effects of the Adler zeros which are important 
in the context of these analyses. Within this parametrization and using ([8]), 

dn 



dE 



{E - 2m w ) 1/2 P a (E)F a (E) = p{E) , 



where, 



P a (E) = AMb 2 /[(M 2 - sY + M 2 T 2 (s)} 



(9) 



(10) 



with s 



E 2 , and 

r 2 (.) : 



s - Am 2 ( s — Sa \ 2 „ , x2 
VE + 2m w (s - s A ) (M 2 
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M 2 



(11) 



1 + 



bi 



M 2 - s A 
b 2 s ( 1 



S) r s - M 



b 2 V M 2 — s s-s A 

Replacing this parametrization of p(E) (with the parameters M, A, sa, b\, b 2 fitted to data 
taken from [20;]) in ([6]), the survival probability is evaluated numerically and is plotted in 
Fig. 1. The different curves in Fig. 1(a) display the contributions of the exponential term in 
()6]), the remaining term (which leads to the power law at large times) and the interference 
of the two terms in the amplitude which appear in P(t). According to 20J|, the a meson 
here has a mass of Er = 542 MeV and a width, = 498 MeV. Though there exist 



+ 



1 

2s 



1 

A 



other predictions 2ll-l23j of the a mass, they all agree on a large width. It is clear from 
the figure that there is a sizable contribution from all terms up to about 15 lifetimes when 
P(t) completely approaches the power law. Fig. 1 (b) makes it clear that the decay law can 
never be approximated by an exponential decay in the case of the a resonance. 

Before proceeding, some remarks regarding the use of such a parametrization are in order. 
Firstly, the calculation of the survival amplitude A(t) requires the analytic continuation of 
the amplitude up to the negative imaginary axis in the lower-right complex energy plane 
(corresponding to the second Riemann sheet in the Mandelstam variable s). This means 
that the knowledge of the amplitude far outside the experimental region is required. It is 
known that there are many parametrizations [22], [23 and references therein] that describe 
the data equally well, but are very different when continued in the complex plane. The 
large uncertainties in the determination of the pole position of the sigma resonance are due 



precisely to this "instability" of analytic continuation. 



detail in 



23]. Secondly, the parametrization of Bugg 20] is valid only up to a region of 



'his issue has been discussed in 



about 1 GeV. As a result of this fact, one encounters several poles in the parametrization at 
high energies which have no physical meaning. Clearly, the occurrence of these poles is an 
artifact of the parametrization and should not be considered in a calculation of the survival 
amplitude of the sigma meson. This is clear alone from the fact that such additional poles 
do not correspond to any known resonant states. Hence, relying on the long energy tail 
of the parametrization, we simply neglect the residues due these poles at high energies. In 
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FIG. 1: Survival probability of the a meson with a width, Tr = 498 MeV as obtained from a 
parametrization of experimental data [^(J. (a) The solid line is the full survival probability P(t), 
the dashed line the contribution of the exponential term, the dashed dotted the remaining part 
and the dotted line is the magnitude of the oscillatory interference term, (b) Comparison of the 
full P(t) with a pure exponential decay law on a log scale. The inlay displays the curves in (a) on 
a linear scale. 

principle, we could have used another parametrization which does not have the drawback 
of such unwanted poles. To clarify this issue in a more detailed way, we refer to the Breit- 
Wigner (B-W) model where we find that the survival probability calculated from the B-W 
model is qualitatively not very different from that obtained using the parametrization in 
[20I (see Fig. 3 to be discussed later). This justifies the neglect of the poles present in the 
parametrization of Ref. [2^ at high energies. 
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FIG. 2: The full survival probability P(t) in a Breit-Wigner model (solid line) as compared to the 
contribution of the exponential term (dashed line), i.e., P(t) = \A exp (t)\ 2 for different values of 
R = Tr/{Er — E t h)- The inlays display the same plots on a linear scale. 

B. Breit- Wigner amplitude 

To get a comparative feeling of the results in Fig. 1 with those of longer lived states, we 
perform some simple model calculations for unstable states with varying lifetimes. We choose 
p(E) to have the standard Breit-Wigner form with a threshold factor and an exponentially 
falling form factor F(E). Thus, 

p B - W (E) — (E — E th y/ 2 x 1/[(E - E R ) 2 + T 2 R / 4] x e~ E / E \ (12) 

where E = 1.1 GeV has been adjusted to match the tail of a realistic parametrization. In 
Fig. 2 we show the plots for unstable states with different ratios R = T R /(E R — E t h) which 
depend on the width as well as the position of the resonance from threshold. It can be seen 
that for narrow states there is a very well defined oscillatory region of transition from the 
exponential to the non-exponential decay law. The oscillatory region shifts to smaller times 
as R increases and for very broad states, the classical approximation of an exponential decay 
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FIG. 3: The full survival probability P(t) in a Breit-Wigner (B-W) model (dashed line) using the 
pole value obtained from Bugg's parametrization is compared in (a) with the result using Bugg's 
parametrization of the amplitude (solid line). In (b) the B-W amplitude is used to compare the 



survival probabilities for different pole values of the sigma meson as given in Ref. 



231 ] by I. Caprini. 



law does not hold good at any time. This is essentially similar to the result shown in Fig. 1 
for the realistic case of the a meson in tttt scattering. Indeed there is no distinct oscillatory 
region of a transition from the exponential to the power law. A similar discussion based on 
a very simplistic model and using a different approach than the one used in this work can 
be found in [24| . 

In Fig. 3a, we compare the full survival probability P(t) evaluated using the Breit-Wigner 



B-W) model (with pole values obtained in 2Qj) with that using the parametrization of Ref. 



20| | . Qualitatively, there is not much difference between the two results, implying that the 
present calculation of the survival probability may not be sensitive to the details of the 
parametrization used. This follows for instance from the fact that the main result of the 
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present work agrees with purely theoretical expectations found in 24|. We therefore have 
to conclude that the overall behaviour of the survival probability as constructed from the 
amplitude is insensitive to the choice of the parametrization. We note that the extra poles 
in Bugg's parametrization are unphysical and an artifact of the parametrization (caused 
probably by the fact that the parametrization is valid only up to a certain energy). The 
comparison with the Breit-Wigner model leads us to the conclusion that we can safely neglect 
these poles which, however, does not imply that both calculations (Breit-Wigner and the 
actual result) are equivalent. That Bugg's model is not a simple Breit Wigner can be seen 
from the expressions. 

Since the B-W model seems reliable once the pole value is given, we can use it to examine 
other parametrizations. Indeed, the survival probability is more sensitive to the pole value 
of the unstable state. This can be seen in Fig. 3b where we use the pole values given in 



231 ] . In 23| I. Caprini performed a detailed analysis of 16 different parametrizations and 



provided certain average best fit pole values for the sigma meson with the correspon ding 
error bars. The various curves presented in Fig. 3b correspond to the pole values from [23 1 
within error bars. It is interesting to note that though the result of a non-exponential decay 
law for the sigma at all times still remains, the behaviour of P(t) in the region where the 
power law sets in, depends on the ratio of the width to the mass of the sigma meson. For 
widths bigger than the resonance mass, the survival probability shows a dip at the onset of 
the power law region. 



IV. SUMMARY 



To conclude, we summarize the findings of the present work: 
[1.] A semi-empirical method to determine the survival probability of an unstable state from 
scattering data is demonstrated with a realistic example of a very broad resonance from 
meson-meson scattering. Though we have found an empirical method of recovering the time 
evolution from experiment, it does depend on the theoretical input of the parametrization 
and the uncertainties associated with it. The power law behaviour at large times which is 
hard to find experimentally is verified for the a resonance. However, the more interesting 
finding is that the decay law for this realistic case of a broad state is never close to an 
exponential. 
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[2.] In order to get a more general view of the behaviour of survival probabilities (P(t)), a 
study using the Breit-Wigner (BW) model was performed and led to the following findings: 

(a) Investigations on the dependence of the critical times (for the transition from an ex- 
ponential to the non-exponential decay) on the positions and widths of the unstable states 
reveal the following: 

(i) for narrow states there exist three distinct regions, namely, an exponential decay, oscil- 
latory transition region and a non-exponential power law. The critical times depend on the 
width as well as the position of the resonance mass from the threshold. They shift to smaller 
values with increasing values of R = T R /(E R — E th ). 

(ii) For very broad states, the decay law does not approach the classical result of an expo- 
nential decay law at any time. Hence, a well-separated exponential followed by an oscillatory 
region does not exist. 

(b) A comparison of P(t) using the BW model and the parametrization in [^(J for the a 
meson shows that overall, P(t) is not sensitive to the details of the parametrization. The 
transition region in P(t) where the power law sets in is sensitive to the ratio of the width to 
the mass of the unstable state. Since the BW result is sensitive only in the transition region 
to the pole values used, the main result that the decay is not exponential at any time still 
remains valid. However, it seems worthwhile to come back in future to the transition region 
which is sensitive to the pole values and use different parametrizations to investigate it. 
[3.] Our method to extract the survival probability has a further significance. The survival 
probability presented in this work is valid only if the system evolving is isolated according 
to its intrinsic dynamics. Interactions with the environment including (repeated) measure- 



ments yield a different picture all-together 



15) . To quote 15j: "The experimentally observed 



survival probability law is exponential at all times. This is due to repeated measurements 
provided At 3> 1 where A is the frequency of the measurements and r the lifetime (for the 



exact definition of the latter see 



151]) . If this is the case, the direct measurement of the sur- 



vival probability defined alone through its intrinsic dynamics can be hampered." However, 
our semi-empirical extraction of this quantity is indirect and does not require a reduction of 
the state. 

These findings should be relevant to most branches of physics where unstable states 
occur. In particular, the example of the a meson presented provides yet another way of 
investigating this elusive scalar meson which has remained a topic of controversy over the 
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years. We note that the a meson, which does not have an exponential decay law at any time 
is an exception among hadron resonances. For all known hadron resonances, the survival 
probability displays an exponential behaviour before the onset of the power law at large 
times. 
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